Phase Controlled Continuously Entangled Two-Photon Laser with 

Double A Scheme 
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We show that an absolute coherent phase of a laser can be used to manipulate the entanglement 
of photon pairs of two-photon laser. Our focus is on the generation of a continuous source of 
entangled photon pairs in the double A (or Raman) scheme. We study the dependence of steady state 
entanglement on the phase and laser parameters. We obtain a relationship between entanglement 
and two-photon correlation. We derive conditions that give steady state entanglement for the 
Raman-EIT scheme and use it to identify region of steady state macroscopic entanglement. No 
entanglement is found for the double resonant Raman scheme for any laser parameters. 

PACS numbers: 



I. INTRODUCTION 



Entangled photon pairs is an integral asset to quantum 
communication technology with continuous variables [l[ . 
A bright source of entangled photon pairs could be useful 
also for quantum lithography . Transient entanglement 
of a large number of photon pairs has been shown to ex- 
ist for cascade scheme an d double Raman scheme 
Q. The transient regime does not provide a continuous 
source of entangled photon pairs that could be as useful 
and practical as typical lasers in c. w. operation. One 
might wonder whether the entanglement still survives in 
the long time limit. 

In this paper, we use the coherent phase of the con- 
trolling lasers to generate a large number (macroscopic) 
of entangled photon pairs in the steady state. This also 
shows the possibility of coherently controlling entangle- 
ment in the steady state. We focus on the Raman-EIT 
scheme (Fig. QJ,) that has been shown to produce non- 
classically correlated photon pairs in single atom and 
many atoms cases. 

First, we discuss the physics of a two-photon emission 
laser using the master equation in Section II. The phys- 
ical significance of each term in the master equation is 
elaborated and related to the quantities of interests (in 
Section III) such as two-photon correlation and Duan's 
Q entanglement measure. In Section IV, we show the 
importance of laser phase for acquiring entanglement. In 
Section V, the steady state solutions for the photon num- 
bers and correlation between photon pairs are given. We 
show that the laser phase provides a useful knob for con- 
trolling entanglement. We then use the results to de- 
rive a condition for entanglement in the double Raman 
scheme. By using proper values of cavity damping and 
laser parameters based on analysis of the condition, we 
obtain macroscopic number of entangled photon pairs in 
the steady state. We also analyze the double resonant 
Raman scheme but find no entanglement. 
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FIG. 1: a) Double Raman atom in doubly-resonant optical 
cavity. The atom is a trapped by an optical dipole force 
and driven by a pump laser and a control laser. Raman-EIT 
(REIT) scheme (fi c , A c (= A) >> fi p ,7) and double resonant 
Raman (DRR) scheme (fi c = Q p , A c = A p = 0) would be the 
focused for analysis, b) Four off-diagonal two-photon emis- 
sion density matrix elements with their respective coefficients 
C k . 



II. PHYSICS OF TWO-PHOTON LASER 

For simplicity, we consider single atom localized in the 
doubly resonant cavity. Using the usual approach we 
obtain the master equation 

■^P = [Ci om i(aipa\ - palax) + C gainl (a\pai - a^p) 
+C\ oss2 (a 2 pa\ - a\a 2 p) + C gain2 (alpa 2 - pa 2 a\) + 



e tlft (Cia 2 pai - C 2 paia 2 + C 3 aipa 2 - C4tii&2/3)]+adj., 

(1) 

with the effective phase iff For double Raman scheme, 
<y?t = <£p + <Pc— (<Ps + fa)- The phases ip a = k a z + cf> a 
of the lasers depend on both the position z of the atom 
and the controllable absolute phases <j> a of the lasers. So, 
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0, the effective 
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FIG. 2: The (Au) 2 + (Av) 2 and g (2) (t) vary in as similar 
manner as a function of time with decoherence 7t, c = 0.67 ac 
(thin line) and without decoherence ybc = (thick line). 
Parameters used are k s = n a = 0.001<7 S , Q. v = 2g 3 ,A p — 
40g 3 ,£l c — 25g s ,A c — 0. We have assumed g a — g s , with 

lac = 7dc = Jab = 7d& = 7- 



phase becomes ipt = P +0 C = 0- The explicit expressions 
for Ciossj, Cgainj and Ck (where j=l,2 and k=l,2,3,4) are 
given in Appendix lAl 

Equation |T]) already includes the cavity damping Li- 
ouvillean hp = — K j{< 1 \< 1 jP + pojctj — 2cijpa ! A. We 

.7=1,2 

find that the relation holds 

C 1 +C 3 = C 2 + C 4 . (2) 

Note that Eq. ([T]) generalizes the master equation for 
the cascade scheme Q in which C3 = C ga i n 2 = 0, and 

Ciossj — foj • 

The Cg a i n j are due to the emissions processes of the 
atom in the excited levels and Raman process via the 
laser fields which provide gain. On the other hand, the 
Ciossj are due to cavity dissipation Kj and absorption pro- 
cesses of the atom in the ground levels which create loss. 
The Ck coefficients correspond to squeezing. Each term 
gives the coherence between rij and rij ± 1 such that the 
difference between the total photon number in the bra 
and in the ket is always 2. Figure [1}d illustrates in two- 
dimensional photon number space the essence of each 
diagonal term in Eq. (fT|). 

The consequence of this relation is that for large 
number of photons rij >> 1, the coherences due to 
a 2 pa\ 1 pa\a 2 ,a\pa 2l h\a 2 p and their adjoint are approxi- 
mately equal. Hence, the contribution of the off-diagonal 
terms vanish and the master equation reduces to the rate 
equation. Since the off-diagonal terms give rise to en- 
tanglement (as we show below), we can understand that 
there will be no entanglement for very large nj. 

III. RELATION BETWEEN ENTANGLEMENT 
AND TWO-PHOTON CORRELATION 

Two-photon correlation for Raman emission doublet 
(RED) (large detuning and weak pump) version of the 
double Raman scheme for single atom [9( and extended 
medium [l(|,[ll| show nonclassical properties such as an- 
tibunching and violation of Cauchy-Schwarz inequality. 



It is useful to show how the nonclassical correlation re- 
lates to entanglement. The normalized two-photon cor- 
relation at zero time delay is 

9^\t) ± + (3) 

- jnwigWit)-!). (4) 

Thus, the g^ 2 \t) does not provide phase 02 1 information 
of the correlation (0.20.1). We introduce the phase via 

(d 2 oi) = I (MO le* 1 . (5) 

Hence the Duan's parameter D((t$.u^ 2 + (Au) 2 can be 
rewritten as 



D(t) = 2[l + ni +n 2 + 2^n 1 n 2 {g ( - 2 K t ) - l)cos0 2 i 

-|<a2)| 2 -|(ai)| 2 -(02)(d 1 )-(a 2 )<« t i>]- (6) 

The terms in the second line would be — |ai| 2 — |a 2 | 2 — 
(aia 2 + a* a 2 ) for an input coherent state. If g^(t) = 
1 we can have entanglement that is independent of the 
phase 02i 

D(t) =2[1- (onaz + ala* 2 )]- (7) 

We shall consider the vacuum state, in which they vanish. 
Clearly, the presence of inseparability or entanglement is 
entirely determined by the phase 02 1 in Eq. ([6]). We now 
find the knob to control entanglement, i.e. cos02i must 
be negative or 7r/2 < 2 i < 37r/2. 

The condition for inseparability or entanglement 
< (Au) 2 + (Av) 2 < 2 can be re-expressed in terms of 
the phase and the two-photon correlation 

ni+n 2 + l ni+n 2 

, < COS 021 < 

2 v / m«2(<? (2) - 1) " 2v/«i« 2 (5 (2) - 1) 

where the lower limit corresponds to maximum entan- 
glement. The midpoint value cos02i = "i+"2+i/2 

gives D(t) = 1. When ni = n 2 . we have — 1 + 1 / 2n < 

cos 02i < 7-^ • Note that for photon antibunching 

9^ ^ 1 the entanglement window for 02 1 can be quite 
large, provided ni,n2 are small too. 

In the limit of large two-photon correlation (pho- 
ton bunching) g^ >> 1 and large photon numbers 
ni,n 2 » 1 the range for entanglement becomes quite 
restrictive. Here, cos 02 1 — , 1 becomes very small 

in magnitude (but negative) and from Eq. ©, D(t) < 2 
the entanglement is small. This explains the results in 
Fig. [2] where large transient correlation is accompanied 
by small entanglement. 

In the long time limit Fig. [5] shows that the corre- 
lation vanishes (corresponds to antibunching) and there 
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is entanglement, D << 2. Thus, both the antibunching 
and entanglement are compatible quantum mechanical 
properties since they manifest in the same way. Duan's 
entanglement increases with time while the correlation 
decreases with time, thus both do not vary in the same 
way. This clearly shows that correlation should be dis- 
cerned from entanglement. However, as expected, the 
decoherence jbc tends to reduce the degree of entangle- 
ment and the magnitude of correlation. 



IV. LASER PHASE FOR ENTANGLEMENT 

Here, we show by using simple example from the reso- 
nant cascde work of Han et. al. @ that the nonzero phase 
of the paired correlation (0,10,2) is necessary for entangle- 
ment. Let us analyze the transient equation (written in 
their notations with zero laser phase) 



d{a\a 2 ) 
dt 



the photon numbers fij do not increase indefinitely. Pa- 
rameters that give non steady state solutions manifest 
as negative value of D and should be disregarded. The 
study of entanglement via nonlinear theory will be pub- 
lished elsewhere. 

In the case of initial vacuum, only 

are sufficient to compute 



dnx dn 2 t^aiag) d(a\a\) 



dt 7 dt 7 dt 



dt 



ia'-a 



Duan's entanglement parameter, where rij 
j = 1,2. From the master equation ((T|), we obtain the 
coupled equations 



^ = n 1 K 1 + e-^(C 1 -C 2 )(o 1 a 2 ) + 



J4i 



(C? - C;)(o\a\) + C gainl + C* inl , (13) 



dn 2 
~~dt 



— — = n 2 K 2 



l0 (C 3 -C 2 )(a 1 a 2 ) 



J4>( 



(C; - C* 2 ){a\a%) + C gain2 + C* (14) 



(9) 



The coefficients for resonant case are such that: fln,(3 22 
are real while (3i 2 — ict\ 2 and /3 2 i = ia 2 i are purely 
imaginary. Clearly we have imaginary value for 

(axa 2 )(t)=i f e- K ^{a 12 (ala 2 )(t') + 



(10) 



a2i((alax)(t') + l)}dt' ~iX 



where K = f3 22 — flu + k 2 + k\ and X is real whose 
expression is not important for the present discussion. 
For initial condition (a,j(Q)) = 0, the Duan's criteria 

D = 2[1 + + (ala 2 ) + (oia 2 ) + (a\al)} (11) 



(a\a 2 ) 



iX-iX} > 2, 



clearly shows D = 2 + 2{(a[a,i) 
there is no entanglement. 

For finite phase <fi associated to the pump laser, the cor- 
relation (&ia 2 ) becomes (aia 2 )e 1 ^ but the photon num- 
bers are not affected. The parameter becomes 



D=2(l 



(a\ai) 



(a\a 2 ) + 2X sin < 



which gives maximum entanglement when 
37r/2, and no entanglement when = 0. 



(12) 
-tt/2 or 



V. STEADY STATE ENTANGLEMENT 

The master equation |T]) is linear and do not include 
saturation. One might wonder whether steady state solu- 
tions. We find that there are steady state solutions when 



(_ _ K 12 )(aia 2 ) = e*[ni(C 3 * - C* 2 )+n 2 {C{ - C* 2 ) - C* 2 ] 



where the gain coefficients/loss are 



(15) 



K ; 



Cgainj + Cgainj — (Clossj + Clossj ) (16) 



K 12 — C gain2 + C* ainl — (C loss2 + C* ossl ). (17) 

The steady state solution for the correlation is (a\a, 2 ) = 
Ee^ where 

E = {-C 32 (K 2 K 1 " 2 + C* 12 C 32 - Ci 2 C3 2 )(Cg ain i + C* ainl ) 
— C* 2 (KxKtf + C 12 C 32 - Ci 2 C 32 )(C gain2 + C* ain2 ) 
+C 2 (K\C\ 2 C% >2 + K 2 C 32 C\ 2 ) — C 2 K\K 2 K'{ 2 



C 2 C* 32 C* 12 {K 1 + K 2 )] 



M 



(18) 



where 



M 



(KxK* 12 + K 2 K 12 ) (C* 12 C 32 ) + c.c. 



- (Ci2C* 32 - C* 12 C 32 ) Z - KxK 2 K 12 Kl 2 . (19) 
The steady state solutions for the photon numbers are 
K 2 K 12 K{ 2 ~ (C\ 2 C 32 K\ 2 + c.c.) 



Hi — (Cgainl + Cgainl)' 

+ (C ga i, l2 + C* in2 )C* 2 Ci 2 



M 

K* 12 + K 12 
M 



+C 2 Ci 2 

+c 2 c* 2 



K 2 K* 12 + (C* 12 C 32 - c.c.) 
M 

K 2 K l2 + {C l2 C* 32 - c.c.) 
M 



(20) 
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n 2 — (Cgain2 + C* in2 ) 



K X K 12 K* 12 - {Cl 2 C 12 Kl 2 + c.c.) g > 2sinfl f and b) if K a n s > S 2 then we need 



+ (Cgainl + Cgainl)C32C32 



M 

K{ 2 + K 12 
M 



+c;q 



2°32" 



KxK* 2 + (C12Q2 - c - c -) 



M 



+c 2 a 



K X K 12 + (C* 12 C 32 - c.c.) 



32" 



M 



(21) 



From Eq. (jlip . the necessary condition for entangle- 
ment is Ee^ + E*e~ l<t> < —(fii + ffa). If £ is real positive 
there would be no entanglement in the region cos <fi > 0. 
Entanglement is still possible even if <j> = provided 
rcal{£?}< 0. Thus, the phase <f> is not necessary for en- 
tanglement, but it provides an extra knob for controlling 
entanglement. 

Let us search for entanglement conditions for the limit- 
ing cases of Raman-EIT (REIT) scheme which produces 
nonclassically correlated photon pairs, and the double 
resonant Raman (DRR) scheme. 



A. Raman-EIT regime 

For this scheme f2 c ,A c (= A) >> il p ,j x (x 
ab, ac, db, dc, be indices for decoherence rates) and A = 0. 
Thus, we have p ba 



in "-{Pbb - Paa) -> 0, p cd = -4^- 



Pdc- It follows from Appendix A that the only finite co- 
efficients in the master equation are 

Clossl — Ks, Closs2 — n a , C 2 — — C4 — iS, (22) 



Ki = -2k s , K 2 = -2 Ka , K 12 = -(k s + K a ) (23) 
where 



- = 9a9s 



A(7 ac7f)C + fi 2 )' 



(24) 



Here, we also have which are used to to write the steady 
state solutions 



n 1 = ~ 2 



n 2 = S 2 



(aia 2 ) =e M HZ 



(k s + K a )[n a n s — z?\ 



(k s + K a )[n a K s ^, 2 ] 

The entanglement criteria can be rewritten as 



(Ks+Ka) 



2sin0, 



< 0. 



(25) 
(26) 
(27) 

(28) 



For negative detuning 5 < 0, there are two possibili- 
ties: a) if K a n s < S 2 entanglement occurs in the re- 
gi° n (k°+k ) 2sin0 t > 5, b) if n a K s > S 2 we have en- 
tanglement in \ 2 sin 6t < 5. Similarly, for pos- 
itive detuning 3 > 0: a) if K a K s < S 2 then we need 



To obtain large entanglement, we tune the cavity 
damping such that \K a K s — C 2 C 2 \ is small and sin t ~ 1. 
Note that the sign of the detuning A in Eq. (f24|) is im- 
portant for entanglement generation. We can arrange the 
signs of these quantities such that Eq. (f28|) is satisfied. 



Figure [3] is plotted using K a = k s = y/l.0l (C 2 C 2 ) and 
A = 407 ac for 6 t ~ 90°. We find that the region tl c ~ A 
gives a large entanglement, but the photon numbers are 
minimum. This prevents having a steady state macro- 
scopic entanglement. We verify that if we change to a 
negative detuning A = — 407 QC there is no entanglement. 

The region of maximum entanglement occur around 
<f> = 90°. Entanglement can occur at a wide range of 
large £l c - However, the photon numbers fij decrease with 
the increase of 51 c . 

Figure d] shows that it is possible to obtain a continu- 
ous bright source of entangled photons. We realize that 
the number of nonclassical photon pairs in REIT scheme 
is limited by the weak pump field. Thus, by increas- 
ing the pump field we can generate more Stokes photons 
(Fig. |4^). At the same time, the detuning is increased 
as well to ensures that the scheme remain in the non- 
classical Raman-EIT regime. By applying the derived 
condition Eq. (|28p we further obtained a larger number 
of entangled photon pairs (Fig. |4}d) . 



B. Resonant Case 



It seems that steady state entanglement in the dou- 
ble resonant Raman case (O c = fi p ,A c = A p = 0) is 
hardly possible. In the following, we investigate this 
analytically. Here, we find C ac ^ c ,C ac . bd ,C bd . a c,C bd . bd 
are real and positive while C ac ,ad, C aCibc , C b d,ad, C bdjbc are 
purely imaginary (positive or negative). Since p cd — 
-i\Pcd\,Pba = ~i\Pba\ all Ci,C 2 ,C 3 ,C4, Ci ossj and C ga inj 
are real and could be negative. So, Kj = 2C ga i n j — 2Ci OSS j 
and K 12 = \{K X +K 2 ). 

For symmetric system, Q p ~ f2 c we find p ab = —p ba = 
Pdc = -Pcd, Pec ^ Pbb and p aa ~p dd [l2j]. Then, C ac . ac = 



bd,bd 



and C a 



,bd 



Q 



If we take T„ 



T, 



dc 



Tdb = 7 (spontaneous decay rate) with T bc = y bc and 
T ad = 27 we further have C bd ,ad = -C aCtad ,C bdibc = 

Cac,bc- 

The steady state solutions for DRR scheme can be 
written as 



rii = n 2 



Cgain(Cgain — Closs) + \C 2 C\ 2 
Cl2 (Cgain Closs )^ 



(29) 



( ai a 2 ) = -e — 2 f— -— (30) 

\ ^12 l^gam ~ Moss; / 

and the entanglement condition with initial vacuum as 
n 1 +n 2 < 2£cos6» t (31) 



a) b) c) 

FIG. 3: Entanglement for Raman- EIT scheme in a) wide view, b) magnified view of focused region (square box) in a), 
and c) mean photon numbers ni (solid line) and n2(dots), where £l p = 7 ac , A p = 407 ac ,A c = 0. Cavity damping values 
K a = k s = a/1.01 (CJC2) ensures a minimum value of the denominator of the condition in Eq. (|28fl . We have assumed 
9a = g s , 76c = with 7 ac = 7 dc = 7 a6 = j db = 7. 




FIG. 4: Macroscopic entanglement for a) larger tt p — 
107a C , A p = 4007ac cavity damping and other parameters the 
same as Fig. [3] b) Q p = 107 ac , A p = 4007 ac but cavity damp- 
ing K a — k s = ^/l.OOl (C2C2). This gives larger mean photon 
numbers and macroscopic entangled photon pairs. 



,-fl 9 2 ^ 2PaQ ~~ Pcc)(fr2p cc + At) - ^-PaaPcc 

(aia2> = e *- 5 

2 7 [^(Pcc-Paa)] 2 -At 2 

(33) 

For strong field, p cc ~ p aa = 0.25. The steady solu- 
tions become fii = n.2 = g^, (0,10,2) = ~£^ e an d 
D = 2(1 + ^7 sin 2 |0t), i.e. no entanglement. 

For weak field, p cc ~ 0.5, p aa ~ 0. The steady solutions 

2 2 ie)^ 2 

are n x = n 2 = f^^r , (0102) = ^j-^zE wi t h K > % 

and hence D = 2(1 + 2 ^ 2g _ g2 sin 2 |^t), again no entangle- 
ment. In the weak field regime, the cavity damping has 
to be sufficiently large to ensure the existence of steady 
state solutions, i.e. At > g 2 /2j. If the cavity damping 
is small At < <? 2 /27 , regions with negative values of D 
and fij would appear which reflect the non-steady state 
regime. 

Thus, we have shown that there is no steady state en- 
tanglement for DRR scheme in both weak and strong 
fields regimes. This is compatible with its classical two- 
photon correlation G( 2 )[II|. However, the REIT photon 
pairs, which are nonclassically correlated, are also entan- 
gled in the steady state. 

VI. CONCLUSION 



where £ is the term in the bracket (...). 

In order to determine whether Eq. (|31[) can be met we 

consider a simpler case where jbc = 0. From Appendix 

2 

B, we have Ci oss - C gain = At and C 12 = ^{pec - Paa), 

2 2 2 

C\ = ^Paa, C 2 = ^(Zpaa ~ Pcc), Goss = (^Pbb + At), 

2 

Cgain — i-Pcc- These results are used to rewrite Eqs. 
(HHJ) and (001) as 

2 

O 2 -p cc K+ ^{2paa- Pcc)(Pcc~ Paa) , x 

"1 = t. t 1 (32) 

2 7 [£{Pcc-Paa)?-K 2 



We have shown that two-photon laser can produce a 
continuous source of entangled photon pairs based on 
the steady state solutions and an entanglement crite- 
ria. We have obtained a relationship between entangle- 
ment and two-photon correlation, and find that both do 
not vary with time in the same manner. We have de- 
rived a condition for steady state entanglement in the 
Raman-EIT (REIT) schemes and showed that steady 
state macroscopic entanglement is possible. We find that 
a large steady state entanglement occurs at the expense 
of smaller photon numbers. We showed that the double 
resonant Raman (DRR) does not generate steady state 
entangled photon pairs for any laser parameters. 
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APPENDIX A: COEFFICIENTS FOR DOUBLE 
RAMAN SCHEME 

The coefficients in Eq. [T]for double Raman scheme are 

Ciossi = \g s \ (Ctd ,adPab + Cbd,bdPbb) + K s, (Al) 

Cgainl = \g S \ 2 {Cbd,bdPdd + CbdMcPdc}, (A2) 

Closs2 — |3a| (Cac.acPcc + C aCl adPcd) + K a , (A3) 
C ga in2 = |ffa| {C ac 

,acPaa + C acM Pba}, (A4) 



Cbd,bc 
Cbd,bd 



I p -Ic 



_?'C) ac a< 
1 ac ± ad ± bc T J P J be T j c J nd 



(A15) 
(A16) 



^ — J ac J ad J 6c J d& "+" 1 P 1 ac 1 ad ' 1 P 1 bc 1 db 

+/ c T a * c T 6 * + 7 c T * d T d6 + (J p - 7 C ) 2 (A17) 



,acPcc + ^acbdPdd + (Cbd,ad + C* cbc )p c a\ A5) 
,acPaa + Cac,bdPdd + C b d,bcPba + C* cbc p(;A6) 
<h = Cbd ,acPaa + ^lc,bdPbb + (Cfed.bc + Cac,ad)Pbc(A7) 
,acPcc 

+ C*c,bdPbb + C*6d .adPcd 

where J^. = <? a , <? s are atom- field coupling strengths, 
C a /3,~fS (a, P,"f,S = a, 6, c, d) are complex coefficients that 
depend on decoherence rates 7 Q/ 3, laser detunings A p , 
A c and Rabi frequencies Cl p , fl c . The p a a,Pab,Pcd (ct = 
a, 6, c, d) are steady state populations and coherences. 



where I p = fi 2 and I c — fl 2 . 
APPENDIX B: COEFFICIENTS FOR RRD 
SCHEME 



From Appendix A, we obtain the coefficients for RRD 
scheme: 

n 2 

C\ = gags-^-2{T bc p cc + T a dPaa} (Bl) 
2 

C 2 = g a g s — ^{T b cPaa + T ad (2p aa -pec)} (B2) 



Cac,ac 
Can. ad 



ad 1 bc 1 db + h 1 ad + 1 c- 1 bc 



•O r bc ^ + Ip-Ic 



(A9) 
(A10) 



12 — ^1 — ^2 — ^3 — ^2 

„2 ^bc + Tad „, 
= ffoffs" ~ 2{p cc -p a 



(B3) 



C a c,bd 



-in. 



- T ad T db + Ip-Ic 



'T 1 * J_ V 1 * 
O O be ~r~ ad 

U - / ■ - u I) 



(All) 
(A12) 



Z = 7 {T ad Tfc c7 + 2/(T a(i + Tb c )}. (B4) 

Taking K s =K s = Kwe also have Ci oss i = Ci oss2 , C gain i = 
C g ain2 and K 2 = Ki = K 12 = 2(C* gain - Ci oss ), so 



Cbd,ac 
Cbd, ad 



= — iQ c 



T bc + T ad 



-T; c r b * c + ip-i c 



(A13) 
(A14) 



2/ITbc , m TT bc 



Closs — \g s \ {—g-Paa + Tad 



Z 



- Pbb ) + K, (B5) 



Cgain — |3s| 2 {Xf 



TT ad + 1 IT ad , , . 

6c 5; Pdd + ^^Pcc } , (B6) 
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